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ABSTRACT 

Using Spitzer Space Telescope and Chandra X-ray Observatory data, we identify YSOs in the Rosette 
Molecular Cloud (RMC). By being able to select cluster members and classify them into YSO types, we 
are able to track the progression of star formation locally within the cluster environments and globally 
within the cloud. We employ nearest neighbor method (NNM) analysis to explore the density structure 
of the clusters and YSO ratio mapping to study age progressions in the cloud. We find a relationship 
between the YSO ratios and extinction which suggests star formation occurs preferentially in the 
densest parts of the cloud and that the column density of gas rapidly decreases as the region evolves. 
This suggests rapid removal of gas may account for the low star formation efficiencies observed in 
molecular clouds. We find that the overall age spread across the RMC is small. Our analysis suggests 
that star formation started throughout the complex around the same time. Age gradients in the cloud 
appear to be localized and any effect the H II region has on the star formation history is secondary 
to that of the primordial collapse of the cloud. 

Subject headings: stars: formation - ISM: clouds - methods: data analysis - individual objects: 
Rosette 



1. INTRODUCTION 

Determining the effect environment has on the star 
forming process is imperative to gaining a full under- 
standing of star formation. Most stars form in embed- 
ded clusters distrib uted throughout a molecular cloud 
(|Lada &; L ada 2003) . These embedded clusters are often 
subject to different environmental conditions based on 
physical their location in the cloud. 

The Rosette Molecular Cloud (RMC) has a ideal lay- 
out for studying star formation and the effects of en- 
vironment. At its northeastern end lies the Rosette 
Nebula which contains the OB association NGC 2244. 
The RMC contains many embedded clusters through- 
ou t the length of t he clo ud. The near-infrared study 
of IPhelps fc Ladal (|1997T) identified 7 embedded clus- 
ters (PL01-07) within the cloud. Subsequent mid- 
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infrared, near-infrared , and X-ray studies have re fined 
and added to this list. Roman-Zuniga ct al. (2008) con- 
firmed the NIR clusters and also discovered two more 
(REFL08-REFL09) using the density distributions of 
near-in frared excess (NIRX) s ource s. Within the PL04 
region, iRoman-Zuniga et al.1 (|2008t ) found the peak of 
NIRX sources ( PL04a) to be spatially coincident with the 
NIR nebulosity. iPoulton et ail (|2008f ) found a concentra- 
tion of Spitzer identified YSOs (cluster D = PL03b) offset 
by 4.3 arcmin (2.0 pc) from the center of the NIRX dis- 
tribution PL03a. Additionally, a distribution of YSOs 
(cluster C = PL02b) was f ound just south of the NIR 
cluster PL02. iWang et all ()2009h found a cluster of X- 
ray sources in the northern end of PL04 indicating the 
presence of a population of Class III sources. The center 
of this distribution (cluster XC = PL04b) is offset from 
the peak of the NIRX distribution by 3.3 arcmin (1.5 
pc). Additionally, a small distribution of X-ray sources 
(XB2) was found spatially coincident with MIR cluster 
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PL02b. 

These embedded clusters are found in different envi- 
ronments. Clusters near the OB association are within 
the ionization front of the H II region. Further away from 
the ionization front, is a region that is associated with a 
high concentration of gas and star formation, where it is 
estimated almo st half of the star formation in the cloud 
is taking place (jRoman- Zufiiga et all [2008). This region 
may have experienced a shock front passing through it. 
Finally at the back of the cloud is a region where two 
embedded clusters have formed beyond the influence of 
t he nebula. 

iRoman-Zuniga et all (|2008[ ) studied the RMC with a 
deep near-infrared JHK survey. In their analysis of NIR 
excess sources they found that the excess fraction within 
the clusters increased with distance away from the the 
center of the Rosette Nebula. This suggested an underly- 
ing age sequence in the cloud where cluster age decreases 
with increasing distance from the nebula. This sequence 
extends beyond the influence of the H II region and it 
was suggested that its origin is from the formation and 
evolution of the cloud. 

In this study we take a closer look at the distributions 
of these young sources and focus on tracking the progres- 
sion of star formation within the cloud. 

2. OBSERVATIONS 

This study makes use of data obtained from the Spitzer 
Space Telescope, the Chandra X-ray Observatory, and the 
FLAMINGOS instrument on the KPNO 2.1 m telescope. 
Figure 1 is a Spitzer IRAC 3-color image of the RMC 
survey region. The blue dashed lines show the extent of 
the FLAMINGOS JHK survey and the red lines show 
the boundaries of the Chandra x-ray observations. The 
locations and names of the embedded clusters are shown. 

2.1. Spitzer IRAC observations and data reduction 

For this study we use IRAC 3.6-8.0 fim and MIPS 24 
/j,m data from program 3394 (PI: Bonnel) available in 
the Spitzer archive. The IRAC mapping covers a total 
projected area of 1.4 x 0.8 degrees 2 for all four IRAC 
bands. Each IRAC pointing has a field of view of 5.12'x 
5.12' with 3 dither positions and an exposure time of 12 
seconds per frame. Additionally, we observed the NIR 
cluster REFL09 with IRAC from a separate program 
40359 (PI: Rieke). The REFL09 IRAC observation (a 
= 06:35:09.11 5 = +03:41:13.7) is a single pointing with 
a field of view of 5.12'x 5.12', 3 dither positions, and an 
exposure time of 12 seconds per frame. 

The IRAC frames were processed using the Spitzer 
Science Center (SSC) IRAC Pipeline vl4.0, and mo- 
saics were created from the basic calibrated data 
(BCD) frames u sing a custom IDL program (see 
IGutermutheFaLl (|2008[ ) for details). The MIPS frames 
were processed usin g the MIPS Data Analysis Tool 
(|Gordon et al.l 12005ft . Source detection and aperture 
photometry was perform using the IDL software pack- 
age PhotVis which is based in part on DAO PHOT 
(jGutermuth et al.l l200i IGutermuth et al.l 12008ft . The 
MIPS 24 /mi frames were processed using the MIPS Data 
Analysis Tool. PSF photometry was performed on the 
MIPS data using the DAOPHOT IRAF package. 

2.2. Chandra observations and data reduction 



The embedded Rosette clusters studied in this paper 
were observed with the Imaging Array of the Chan- 
dra Advanced CCD Imaging Spectrometer (ACIS-I; 
iGarmire et al.l (|2003[ )) onboard of the Chandra X-ray Ob- 
servatory, which has a 17' x 17' field of view in a single 
pointing. These observations were taken on 2010 Decem- 
ber 3 (ObsID 12388, covering the PL3 cluster), 2010 De- 
cember 10 (ObsID 12142, PL6), 2011 January 14 (ObsID 
12387, PL1), and 2011 January 18 (ObsID 12386, REFL 
9), with net exposure time of 19.6 ks, 39.3 ks, 24.5 ks, and 
34.6 ks, respectively. They significantly supplement the 
previous Chandra campaign of the Rosette comp l ex re- 
porte d inlTownslev et all (|2003) and iWang eFaLl (|2008l 
2009,120100 which consisted of four ~20 ks ACIS-I snap- 
shots in January 2001, a deep 75 ks ACIS-I image in 
January 2004 centered on the 05 star HD 46150 in NGC 
2244, and one 20 ks ACIS-I pointing at the NGC 2237 
sub-cluster in 2007. All images were taken in "Timed 
Event, Very Faint" mode (5 pixel x 5 pixel event is- 
lands) . 

We follow the same c u stomiz e d da ta reduction 
described in IWang et al.l (|2007l I2008D using the 
Chandra Interactive Analysis of Observations (CIAO, 
iFruscione etaLl ((2006); version 4.3) package provided by 
the Chandra X-ray Center. A detailed description of the 
X-ray data analysis will be presented in a separate pa- 
per (Wang et al. 2013, in preparation). For each of the 
ACIS fields, we extracted X-ray images in the 0.5-7 keV 
band, and applied the sour ce detection algorithm wavde- 
ccman et~aT]|2002]) with a range of wavelet scales 

(from 1 to 16 pixels in steps of \/2) and a source signif- 
icance threshold of 1 x 10" 5 to produce a list of candi- 
date sources. X-ray event extraction was m ade with our 
custo mized IDL script ACIS Extract (AE; iBroos et al.1 
2010). Using the Miscalculated probability P B that 
the extracted events are solely due to Poisson fluctua- 
tions in the local background, we rejected sources with 
Pb > 0.01, i.e. those with a 1% or higher likelihood of 
being a background fluctuation. The trimmed source list 
includes 431 valid X-ray sources. 

2.3. Near-infrared photometry data 

This study makes use of photometry data from 
the FLAMINGOS JHK im aging survey of the RMC 
(jRoman-Zuniga et al.ll2008ft which is part of the NO AO 
survey program "Toward a Complete Near-Infrared 
Spectroscopic and Imaging Survey of Giant Molecular 
Clouds" (PI: E. A. Lada). The completeness limits for 
the survey are J = 17.25, H = 18.00, and K = 18.50. 

3. ANALYSIS 

3.1. Dust extinction distribution in the RMC 

A dust e xtinction map w as made with the NICEST 
algorithm (|Lombardil 12009ft on the FLAMINGOS pho- 
tometry catalog. We removed Class 0/1 and Class II 
sources from the catalog as those sources have large in- 
trinsic red colors and tend to bias the map near the 
centers of clusters. The extinction is estimated toward 
background sources by comparing the source NIR col- 
ors with the intrinsic color distribution measured from a 

2 http: //cxc. harvard. edu/ciao/threads/wavdetect/ 

li http : //www . astro .psu. edu/xray/docs/TARA/ae_users_guide .html 
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Figure 1. Spitzer 3-color image of the Rose tte Molecular Cloud [4.5 fim (blue), 5.8 /an (green), 24 /an (red)]. Blue dashed lines show the 
boundaries of the FLAMINGOS JHK survey IjRo man-Zufiiga ct al. 2008). Red lines show the boundaries of the Chandra x-ray observations 
HWang et al.H200g . this paper). 



nearby extinction-free control region. The control field 
is located at a similar galactic latitude as the RMC and 
was selected from an IRAS 2 5 um map as being devoid 
of dust thermal emission (see lRoman-Zuniga et al.ll2008l 
Fig. 2). Spatial smoothing is then applied to the extinc- 
tion values to create the extinction map. The smoothing 
creates a map of the weighted mean as a function of po- 
sition on the sky. The weighting function is composed of 
a Gaussian which weighs the individual extinction mea- 
surements based on angular distance from the center of 
the map point plus a correction term which compen- 
sates for the bias due to the low numbers of background 
sources at high extinctions. In the case of the RMC, 
the optimal Gaussian width was found to be 90". This 
is capable of resolving column density structures of the 
molecular cloud with projected sizes of about 0.5-1.0 pc. 

Figure 2 shows the extinction map of the RMC. The 
contour levels indicate Ay = 8, 10, 12, 14, 16, 20 mag. 



The locations of the embedded clusters are indicated. 
Comparison of the location of the embedded clusters with 
the extinction distribution suggest a correspondence be- 
tween star formation with the highest extinction regions 
of the cloud. 

3.2. Identification of YSOs 

Spitzer and Chandra ACIS data were used to identify 
Young Stellar Objects (YSOs) in the Rosette Molecu- 
lar Cloud. In order to select the Cl ass 1/0 and Class 
II sou rces, we employ the color cu ts of iGutermuth et all 
(2008) and lKrvukova et all (|2012l ) to our Spitzer catalog. 
For objects which have IRAC (3.6 ^m, 4.5 ^m, 5.8 fim, 
8.0 /im) and MIPS 24 /im detections, we identify Class 
1/0 sources using the following criteria: 

[4.5] - [5.8] > 1 
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Figure 2. NICEST extinction map of the Rosette Molecu- 
lar Cloud. The contours represent extinction levels Ay = 
8,10,12,14,16,20 mag. Location of embedded clusters are indicated. 

[4.5] - [5.8] > 0.7 and [3.6] - [4.5] > 0.7 

and 

[4.5] - [24] > 4.7 

From the remaining catalog we select Class II sources 
that have all of the following criteria: 

[4.5] - [8.0] > 0.5 
[3.6] - [5.8] > 0.35 
[3.6] - [5.8] < 3.5Q4.5] - [8.0]) - 1.75 

The Gutermuth [4.5]-[5.8] color cut that distinguishes 
between Class 1/0 and Class II sources is particularly 
useful as the extinction curve is relatively flat in that 
wavelength range and thus the color cut is insensitive to 
extinction. For sources without an IRAC 4.5 /xm detec- 
tion, we select Class 1/0 sources by 

[5.8] - [24] > 4.5 and [24] < 6 

and for Class II sources 

[3.6] - [5.8] < 0.35 and 2.0 < [5.8] - [24] < 4.5 

For objects that do not have a MIPS 24 /im detection, 
we use the previous IRAC color cuts for objects with a 
[4.5]-[5.8] color and the additional requirement of 

[5.8] - [8.0] < 1 

in order to filter out AGN and PAH galaxies. 

Class III sources do not display significant infrared ex- 
cess and thus need to be identified another way. Fortu- 
nately, YSOs are known to emit X-rays at levels that can 
range many orders of magnitude above main s equence 
stars (jPreibisch et alj|2005t iFeigelson et al.ll2007t) . Thus, 
X-ray observations can efficiently identify YSOs in molec- 
ular clouds. Class III candidate sources were selected 
from our Chandra ACIS ob servations and the previously 
published X-ray catalog of W ang et al.l (|2009f ). Sources 
with colors of Class 1/0 or Class II objects were then re- 
moved to create a catalog of Class III sources. In order 
to deal with extragalactic contaminants in our sources 
list, we also removed sources that did not have a NIR 



counterpart in the FLAMINGOS catalog. In our anal- 
ysis of the X-ray properties, we found that the sources 
without a NIR counterpart had an average hardness ra- 
tio consistent with extragalactic sources (Wang et al., in 
preparation) . 

3.3. Spatial distribution of the different YSO classes 

By analyzing the distributions of YSO classes sepa- 
rately one can probe the evolution of star forming re- 
gions. The different YSO classes represent different evo- 
lutionary stages, with Class and I sources representing 
the youngest sources still embedded in their envelopes, 
Class II sources are a later stage of sources still accret- 
ing material from their disks, and Class III sources are 
diskless pre-main sequence stars. 

We employed the k-nearest neighbor (k-NN) density 
estimation algorithm, often referred to as the nearest 
neighbor method (NNM), to analyze the structures and 
distributions of the YSO classes. This method allows us 
to map the density distri butions and subsequently iden- 
tify regions of clustering (jCasertano fc Hutl lT985). The 
algorithm measures the distance, Dk, between a source 
point and its k th nearest neighbor and estimates the lo- 
cal stellar number density as Hk = (k — V)lnD\. For 
the Class 1/0 objects we used a k=5, and for the Class II 
and Class III sources we used k=10. Using k values larger 
than one has the advantage of lessening the influence pos- 
sible non-members of the set, e.g. background galaxies, 
have on the local density estimate (|Ferreiral2010|) . This is 
important when using Spitzer data as many background 
galaxies have colors that are simi lar to those of YS Os in 
both the near- and mid-infrared (F oster et aLl feOOS). 

Figure 3 presents the NNM density maps of the Class 
1/0, Class II, and Class III sources in the RMC. The 
dashed lines show the spatial boundaries of the survey 
data used to identify the sources. The embedded clus- 
ters are spatially coincident with the high stellar density 
regions in the maps. Using the NNM maps of the Class 
II and Class III sources, we have identified a new cluster 
(PL06b). The maps reveal that the different YSO classes 
can have very different spatial distributions. For exam- 
ple the Class 1/0 and Class III concentrations appear to 
be separate while the Class II sources appear more wide 
spread. This is especially evident in the central region of 
the cloud. Figure 4 shows the NNM density maps of the 
central region for the Class 1/0 and Class III sources. A 
high density of Class 1/0 sources is found in the center 
of this region and is spatially coincident with the highest 
extinction. In contrast the older Class III sources are 
found primarily around the perimeter on north and east 
sides of the central region. Class II sources are found 
throughout the central region. By numbers, the Class II 
sources dominate the cloud. The different density distri- 
bution of the YSO classes in the central region suggest 
a progression of star formation from the outer north and 
east sections towards the central. 

3.4. Age distribution through YSO ratios 

In order to study the progression of star formation in 
the RMC, we investigate how the ratios between differ- 
ent YSO classes change throughout the cloud. As time 
progresses, Class 1/0 sources will lose their envelopes and 
evolve into Class II sources. Because the lifetime of the 
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Figure 3. Nearest-Neighbor (NNM) density maps of the RMC 
for Class 1/0, Class II, and Class III sources. The density contours 
are fj, = 2.6,4.6, 7.7, 12.8 stars pc — 2 . The background contours are 
Ay = 8, 10, 12, 14, 16, 20. Dotted lines in the first two panels show 
the coverage of the Spitzer observations and the dotted lines in the 
last panel show the coverage of the Chandra observations. 

Class I/O phase is less than the Class II phase, the ratio 
of Class I/O to Class II sources s hould decrease with the 
age of the region fcf. lMversll2012l ). Similarly, as time pro- 
gresses, Class II sources will lose their disks and evolve 
into Class III sources. Disk fracti ons in clusters hav e 
an empirical relationship with age (|Haisch et al.l f2001). 
From different Spitzer IRAC studies of young clusters, 
the average disk fractions at different cluster ages are 
75% at 1 Myr, 50% at 2-3 Myr, 20% at 5 Myr, and 5% 
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Figure 4. Nearest-Neighbor (NNM) density maps of the RMC 
main core region for Class 1/0 and Class III sources. The density 
contours are fj, = 2.6,4.6,7.7,12.8 stars pc -2 . The background 
contours are Ay = 8, 10, 12, 14, 16, 20 



Table 1 



Age (Myr) 


D.F. 


Rii-.iii 


1 


75% 


3.0 


2-3 


50% 


1.0 


5 


20% 


0.3 



Note. — Columns 1 and 2 are cluster ages and average disk 
fractions from the literature (Williams & Cicza 2011). Column 3 
is the associated Class II to Class III ratio. 



at 10 Myr ()Williams fc Ciezal |20TTI ) . Thus, the Class II 
to Class III ratio within a region will also decrease with 
age (Table 1). This relationship between the Class II to 
Class III ratio and age allows us to probe age progres- 
sions within the cloud. 

We constructed ratio maps by estimating the ratio of 
the number of one class of sources, N\, to that of another, 
N 2 , (e.g. Class 1/0 to Class II) within a projected region 
of the sky across a grid . We chose the size of the region 
by trying to minimize both resolution and uncertainty. 
The probability density function of the ratio, R, is 



p(R\N u N 2 ) = 



R N ^(N 1 +N 2 + 1)1 



(R+ l)Nx+N 2 +2 jvi!JV 2 ! 

assuming a uniform prior for the ratio (jJin et al.ll2006h . 
We use the expectation value of the ratio as our estima- 
tor, 

JVi + 1 



R 



N 2 



with variance, 



4 



(jVi + l)(7Vi + JV 2 + l) 
Ni(N 2 -l) 



In each point on the grid, for a non-zero ratio value, we 
require aji/R < 0.50. For the Class 1/0 to Class II ratio 
we chose a circular region with radius 1.1 pc and for the 
Class II to Class III ratio a circular region with radius 
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1.4 pc. The radius for the Class II to Class III ratio is 
larger because the magnitude cuts (see next paragraph) 
reduce the number of sources. The regions are sampled 
across the grid at intervals of one-third the region radius. 
The created ratio map provides a visual representation 
of age gradients and can then map the progression of star 
formation. 

In order to compare Class II and Class III sources it 
is necessa ry to make sure the sam ples are uniform and 
unbiased (jGutermuth et all I2004D . We first restrict the 
samples to sources which are detected in the FLAMIN- 
GOS JHK survey. We estimate the Class II sample is 
complete to H = 15 and the Class III sample to be 
complete to H = 14. We use the H-band luminosity 
as the emission from disks at this wavelength is negli- 
gible. We de-redden our samples to a 2 Myr isochrone 
(jBaraffe et al.l 119981 ) and limit the samples to H < 14, 
making the samples complete to the same depth across 
the cloud. 

Figure 5 shows the two ratio maps, the Class II to 
Class III on the left and the Class I/O to Class II on 
the right. These maps are over-plotted on top of the 
extinction contours. For the Class II to Class III ratio, 
the contour levels are Rh-jii = 0.5, 1.0, 2.0, 3.0. For the 
Class I/O to Class II ratio, the contour levels are Riji = 
0.1,0.3,0.4,0.5. Visual inspection of the maps suggest 
a correlation between the ratios and extinction, where 
higher ratio values are spatially correlated with higher 
extinction. 

To further investigate the correlation between these ra- 
tios and extinction we also determine the average extinc- 
tion in each region. We binned the regions by extinction 
into 1 mag bins and calculated the weighted mean ra- 
tio in each bin. Figure 6 presents a plot of extinction 
bin versus mean ratio of Class II to Class III sources. 
The figure shows that the ratio increases monotonically 
with extinction. The plot can be fitted with a shallow 
linear fit for Ay < 13 mag, and a steeper linear fit for 
Ay > 13 mag. Thus, decreasing age is related to in- 
creasing extinction. Figure 7 shows the extinction bin 
versus mean Class 1/0 to Class II ratio. At high extinc- 
tions, Ay > 17 mag, the Class 1/0 to Class II ratio also 
steeply increases with extinction. However, this ratio ap- 
pears flat for the extinction range Ay = 8-17 mag. To 
investigate this further we look at the contributions to 
the graph from the main core and cluster PL07 sepa- 
rately (Fig. 7 b & c). We find that for the main core 
of the cloud, where most of the star formation is taking 
place, the ratio has a positive monotonic relation with 
extinction. Cluster PL07 region is different; It appears 
to have a relatively flat relationship between ratio and 
extinction, with a possible increase at Ay = 7-8 mag. 
There is a small offset between the extinction peak and 
Class 1/0 density peak; which suggests a recent expul- 
sion of gas from the center of PL07. For the cloud as a 
whole, the mean ratios between the Class 1/0 and Class 
II sources have a small range (0.20-0.60) which may cor- 
respond to a narrow range of ages or perhaps evidence 
for continuous star formation. 

The relationships between the two ratios and extinc- 
tion reveal a trend of decreasing age with increasing ex- 
tinction. This implies that stars may be preferentially 
forming in the highest extinction parts of the molecular 
cloud. 



4. DISCUSSION 

4.1. Star formation and column density 

The Class II to Class III ratio map traces regions with 
estimated ages of 0.5-3 Myr, and these regions are spa- 
tially coincident with extinction values of Ay = 4-18 
mag. Using the empirical relationship between disk frac- 
tion and cluster age, we can fit a power law relation to the 
Class II to Class III ratio and age. Then we can directly 
study the relation between age and extinction. Figure 8 
shows a plot of estimated age versus extinction. The plot 
suggests that the column density of gas decrease expo- 
nentially with time above Ay = 5 mag, with a half-life, 
ii/2 = 0-4 Myr. We estimate the average rate at which 
the column density of gas decreases as 

S - 10~ 4 M yr" 1 pc" 2 

which is over an order of magnitude larger than the 
star formation rate measured in nearby molecular clouds 
(|Lada fc Ladal [20031: lEvans et al J 12009 ). Thus most of 
this gas is not removed through formation of stars and 
is possibly being relocated to other regions of the cloud. 
This is co nsistent with the study of nearby molecular 
clouds by lLada et "all (|2010[) which demonstrated that 
star formation on the scale of few Myr has a negligible 
effect on the total mass of the cloud. 

In the RMC main core, we find a similar relationship 
between the Class 1/0 to Class II ratio and extinction. 
This suggests that star formation occurs preferentially 
in high extinction regions. Additionally, we find that 
over half of the clustered (n > 1.3 stars pc -2 ) Class 1/0 
sources are found at Ay > 15 mag and all of them at 
Ay > 7.5 mag. This too i s consistent with t he study of 
nearby molecular clouds bv lLada et al.1 (|2010l ) that found 
an extinction threshold of Ay ~ 7 mag, above which the 
star formation rate was proportional to the mass of the 
cloud measured above that threshold. 

4.2. Cluster Properties 

The relationship between star formation and column 
density may have consequences for the formation and 
evolution of clusters. Table 2 shows the properties and 
YSO content of the clusters. The stellar content is mea- 
sured within a 1 parsec (2.14') radius of the listed cluster 
center. 

Most of the embedded clusters in the RMC have a 
Class II to Class III ratio, Rnju, between 1 and 2. 
The weighted mean ratio of all the embedded clusters is 
Riijii = 1.2 ± 0.2, suggesting that most of the clusters 
started forming around the same time. 

There is a group of clusters with ratios suggesting a 
more recent episode of star formation, having Rum > 
3.0. This group includes clusters PL2a, PL02b, PL03b, 
and REFL09. These clusters however do not appear to 
have a significant Class 1/0 population. This suggests 
that the star formation in these clusters is more coeval. 

The cluster REFL08 appears to have the most recent 
episode of star formation. This cluster has Class 1/0 
sources concentrated in dense filamentary structures seen 
in extinction. It has the highest ratio of Class II to Class 
III sources which is consistent with an age less than 1 
Myr. Our X-ray observations do not cover embedded 
cluster PL07. This cluster has the second highest con- 
centration of young Class 1/0 objects. However, without 
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Figure 5. Left: Class II to Class III ratio map. The contour levels are Rzi-.Iir = 0.5, 1.0,2.0,3.0. Right: Class 1/0 to Class II ratio map. 
The contour levels are Rj.jj = 0.1,0.3,0.4,0.5,0.7. The background contours are Ay = 8,10,12,14,16. Dotted lines show the coverage of 
the Spitzer observations and dash dotted lines show the coverage of the Chandra observations 
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Figure 6. Plot of extinction bin versus mean Class II to Class 
III ratio. 



knowledge of its Class III content, an age estimate for 
this region is not possible. 

Although there are some differences between the clus- 
ters, the age spread of these clusters is nonetheless small. 
The ages inferred from the YSO ratios for most of the 
clusters are between 1 to 3 Myrs. 

4.3. Age gradients in the RMC main core 

The RMC main core is composed of clusters PL04(a & 
b), PL05, and REFL08. Each of the clusters is associ ated 
with a CO clump identified in lWilliams et al.l (|1995[ ) CO 
survey of the RMC (Table 3). These three clusters are 
characterized by having more YSOs than the other clus- 
ters in the cloud, wh i ch is consistent with the study by 
IRoman-Zuniga et al.l (|2008) where it is estimated that 
half of the star formation in the whole cloud happens in 
this region. 

The REFL08 sub- region has the highest density of pro- 
tostars and is spatially coincident with the gas surface 
density peak. Its ratio of Class 1/0 to Class II sources 



o 

CO 



CO 

JS 
O 



CO 

JS 
O 



a) whole cloud 



I 



~i 1 1 1 1 1 r 

b) main core 



I 



i 1 1 1 1 1 r 

c) PL07 



I 



mi 



Figure 7. Plot of extinction bin versus mean Class 1/0 to Class 
II ratio. 
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Table 2 

Cluster Properties 



Cluster 


R.A. 


Dec. 


I/O 


II 


III 


II 

H < 14 


III 

H < 14 


Rii-.iii 


Ay 


PL01 


97.96 


4.32 


5 


11 


11 


10 


8 


1.4±0.7 


8.1±1.4 


PL02 


98.31 


4.59 


2 


10 


6 


7 


2 


4.0±4.0 


6.4±2.2 


PL02b 


98.31 


4.53 


3 


15 


5 


10 


3 


3.7±2.9 


7.7±1.0 


PL03a 


98.38 


4.00 


6 


■1 


8 


2 


6 


0.5±0.4 


10.5±2.5 


PL03b 


98.42 


4.06 


2 


19 


5 


12 


1 


3.3±2.2 


9.9±1.9 


PL04a 


98.53 


4.42 


7 


25 


14 


18 


10 


1.9±0.8 


10.7±2.4 


PL04b 


98.55 


4.47 


2 


32 


30 


21 


17 


1.3±0.4 


7.2±0.8 


REFL08 


98.56 


4.32 


18 


30 


4 


19 


1 




17.8±2.4 


PL05 


98.63 


4.32 


2 


32 


15 


22 


11 


2.1±0.8 


11.3±2.2 


PL06 


98.66 


4.21 


4 


10 


13 


9 


9 


1.1±0.5 


10.8±1.6 


PL06b 


98.59 


4.20 





12 


10 


7 


7 


1.1±0.7 


10.2±1.2 


REFL09 


98.78 


3.71 


5 


36 


4 


16 


3 


5.7±4.4 


16.6±2.7 


PL07 


98.88 


3.98 


12 


32 


ND 








10.6±1.4 



Note. — The number of sources are counted within 1 pc (2.14') of the cluster centers 



Table 3 

Associated CO clumps 



I 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Estimated Age (Myr) 

Figure 8. Plot of estimated age versus mean extinction. The age 
is estimated from the Class II to Class III ratio. The red dotted 
line shows the exponential fit with half-life, t\i% = 0.4 My, above 
Av = 5 mag. 



and its dearth of Class III sour ces suggest the age of this 
region to be less than 1 Myr. iHennemann et afl ()2010f ) 
using Herschel observations found 27 protostars in this 
region, 7 of which were classified as very young Class 
candidates. This region appears filamentary in the MIR 
Spitzer images with the main filament running north to 
south. The southern end of the main filament is coinci- 
dent with the center of NIR cluster REFL08 and is at 
the intersection of two smaller filaments. This region ap- 
pears to be the youngest region of the cloud. Cluster 
PL04 is located north of REFL08, while PL05 is located 
to the east. Both clusters, PL04 and PL05, appear to 
be older with Class II to Class III ratios consistent with 
ages 1-2 Myrs. 

The age gradient across this region as a whole suggests 
that star formation began in clusters PL04 and PL05 
first, followed by star formation in REFL08 as a more 
recent event. Based on the YSO content, we estimate 
the age difference across the region to be about 1 Myr. 
It is possible that star formation feedback from clusters 
PL04a and PL05 pushed the gas into its present state 
and thus triggered the formation of REFL08. 

Alternatively, the formation of cluster REFL08 may 
not have been triggered by PL04 and PL05. The age 



Cluster 


Clump 


V 


Av 


PL04 


1 


15.6 


2.1 


PL05 


19 


12.2 


1.7 


REFL08 


6 


10.8 


1.9 



Note. — CO clump data from Williams et al. (1995) 

progression may be a consequence of the formation of 
the cloud. This scenario is consistent with the numeri- 
cal simulations o f dyna mic molecular cloud formation by 
lHartmann et al.l ([2012). In these simulations, stars can 
form from the initial density fluctuations resulting from 
turbulence during the formation of a molecular cloud. 
These stars form before the global gravitational collapse 
of the cloud leads to a main phase of star formation. 

4.4. Star Formation as a function of distance from 
NGC 22U 

We find that the spatial distribution of young Class 1/0 
sources throughout the cloud has little or no correlation 
to location in the cloud or distance form NGC 2244. The 
age of NGC 2244 is estimated to be 2-3 Myrs which is 
consistent with the ages of the other clusters inferred 
from the YSO ratios. 

Age progressions appear to be limited to small regions, 
where surface density enhancements lead to more recent 
star formation episodes. Clusters PL02a and PL02b, spa- 
tially coincident with the visible edge of the H II region, 
have Class II to Class III ratios consistent with a re- 
cent (~ 1 Myr) episode of star formation. The forma- 
tion of these clusters may have been trigge red by NGC 
2244, consistent with the X-ray analysis of iWang et al.1 
(2009). Our analysis finds a local age progression sur- 
rounding REFL08 which is consistent with the studies of 
iRoman-Zuniga et al.l ([2008) . This progression does not 
extend through the cloud. It appears that any effect the 
H II region has on the star formation history is secondary 
to that of the primordial collapse of the cloud. Our anal- 
ysis suggests that star formation started throughout the 
complex due to a primordial collapse and has progressed 
through a series of localized episodes of formation closely 
following each other. 
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5. CONCLUSION 

In our analysis of the different YSO classes in the RMC 
we find that they often have different density distribu- 
tions. We use the YSO ratios to study age gradients 
across the cloud to better understand how star forma- 
tion has progressed. The relationships between the YSO 
ratios and extinction suggest that star formation in the 
cloud occurs preferentially in high extinction regions and 
that the column density of gas rapidly decreases as the 
region evolves. This suggests rapid removal of gas may 
account for the low star formation efficiencies observed 
in molecular clouds. We find that progressions of star 
formation appear to be localized with a small overall age 
spread across the cloud, consistent with star formation 
starting across the cloud at roughly the same time. The 
distribution of YSOs in the RMC show little or no effect 
of NGC 2244 on the relative ages of the clusters except 
for the possible triggering of clusters PL02a and PL02b. 



We thank Charles Lada for useful discussions and sug- 
gestions in the writing of this manuscript. We also thank 
Rob Gutermuth for help with the Spitzer IRAC photom- 
etry. This work is based in part on archival data ob- 
tained with the Spitzer Space Telescope, which is oper- 
ated by the Jet Propulsion Laboratory, California Insti- 
tute of Technology under a contract with NASA. Sup- 
port for this work was provided by an award issued 
by JPL/Caltech, a NASA LTSA Grant NNG05GD66G, 
and a NSF grant AST-II09679. Partial support for this 
study provided by the Chandra ACIS Team through the 
SAO grant SV4-740I8. J.Y. thanks his GSRP advisor, 
Malcolm Niedner, for encouragement and support. J.Y. 
also acknowledges support from a NASA Goddard Space 
Flight Center GSRP fellowship NNX10AM07H and a 
Space Grant Fellowship from the Florida Space Grant 
Consortium. C.R.-Z. acknowledges support from CONA- 
CyT, Mexico, project I52I60 and from project UNAM- 
DGAPA-PAPIIT IA101812. J.W. acknowledges support 
from NASA grant GO1-12030X. 

Facilities: Spitzer (IRAC,MIPS) CXO (ACIS) . 

REFERENCES 



Baraffc, I., Chabrier, G., Allard, F., et al. 1998, A&A, 337, 403 
Broos, P. S., Townsley, L. K., Feigelson, E. D., et al. 2010, ApJ, 
714, 1582 

Casertano, S., & Hut, P. 1985, ApJ, 298, 80 
Evans, N. J., II, Dunham, M. M., j0rgensen, J. K., ct al. 2009, 
ApJS, 181, 321 

Feigelson, E., Townsley, L., Giidel, M., et al. 2007, Protostars and 

Planets V, 313 
Ferreira, B. 2010, Ph. D. thesis, Univ. of Florida 
Foster, J. B., Roman-Zuniga, C. C, Goodman, A. A., ct al. 2008, 

ApJ, 674, 831 
Freeman, P. E. ct al. 2002, ApJS, 138, 185 
Gordon, K. D., Rieke, G. H., Engelbracht, C. W., et al. 2005, 

PASP, 117, 503 
Fruscione, A., et al. 2006, Proc. SPIE, 6270 
Garmire, G. P., Bautz, M. W., Ford, P. G., ct al. 2003, 

Proc. SPIE, 4851, 28 
Gutermuth, R. A., Megeath, S. T., Muzerolle, J., et al. 2004, 

ApJS, 154, 374 
Gutermuth, R. A., et al. 2008, ApJ, 674, 336 
Haisch, K. E., Jr., Lada, E. A., & Lada, C. J. 2001, ApJ, 553, 

L153 

Hartmann, L., Ballesteros-Paredes, J., & Heitsch, F. 2012, 

MNRAS, 420, 1457 
Hennemann, M., Motte, F., Bontemps, S., et al. 2010, A&A, 518, 

L84 

Jin, Y. K., Zhang, S. N., & Wu, J. F. 2006, ApJ, 653, 1566 
Kryukova, E., Megeath, S. T., Gutermuth, R. A., ct al. 2012, AJ, 
144, 31 

Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57 

Lada, C. J., Lombardi, M., & Alves, J. F. 2010, ApJ, 724, 687 

Lombardi, M. 2009, A&A, 493, 735 

Muench, A. A., Lada, C. J., Luhman, K. L., et al. 2007, AJ, 134, 
411 

Myers, P. C. 2012, ApJ, 752, 9 

Phelps, R. L., & Lada, E. A. 1997, ApJ, 477, 176 

Poulton, C. J., Robitaille, T. P., Greaves, J. S., ct al. 2008, 

MNRAS, 384, 1249 
Preibisch, T., Kim, Y.-C, Favata, F., et al. 2005, ApJS, 160, 401 
Roman-Zuniga, C. G., Elston, R., Ferreira, B., ct al. 2008, ApJ, 

672, 861 

Townsley, L. K., Feigelson, E. D., Montmerle, T., et al. 2003, 
ApJ, 593, 874 

Wang, J., Townsley, L. K., Feigelson, E. D., ct al. 2007, ApJS, 
168, 100 

Wang, J., Townsley, L. K., Feigelson, E. D., ct al. 2008, ApJ, 675, 
464 

Wang, J., Feigelson, E. D., Townsley, L. K., et al. 2009, ApJ, 696, 
47 

Wang, J., Feigelson, E. D., Townsley, L. K., et al. 2010, ApJ, 716, 
474 

Williams, J. P., Blitz, L., & Stark, A. A. 1995, ApJ, 451, 252 
Williams, J. P., & Cieza, L. A. 2011, ARA&A, 49, 67 



